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Abstract 

We present evidence for the universality of the continuum limit of the 
scale dependence of the renormalization constant associated with the op- 
erator corresponding to the average momentum of non-singlet parton den- 
sities. The evidence is provided by a non-perturbative computation in 
quenched lattice QCD using the Schrodinger Functional scheme. In par- 
ticular, we show that the continuum limit is independent of the form of 
the fermion action used, i.e. the Wilson action and the non-perturbatively 
improved clover action. 
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We presented in ref. |L] a quenched non-perturbative calculation for the renor- 
malization constant associated with the operator corresponding to the average 
momentum of non-singlet parton densities. The calculation was based on a fi- 
nite size recursive scheme that allows us to reconstruct in the continuum the 
non-perturbative evolution of the renormalization constant in the Schrodinger 
Functional (SF) scheme @, |3|, [|, [|. In ref. |IJ we used a standard Wilson action 
where the continuum limit is approached with a rate that is linear in the lattice 
spacing a in leading order. The results for the step scaling functions presented in 
ref. |TJ also showed sizeable quadratic corrections that had to be included in the 
fit and made the continuum extrapolation more hazardous. 
We repeated the same calculation with a non-perturbatively improved clover ac- 
tion [|J that changes the form of the lattice artefacts. Linear corrections are still 
to be expected, because the 0(a)-improvement of the action should be comple- 
mented with the improvement of the operators and of the boundary counterterms 
in order to lead to a full cancellation of effects appearing linear in a. However, 
the continuum limit cannot depend upon the discretization chosen and should be 
universal. Our results do indeed show such a universality and therefore put our 
previous continuum extrapolations on a firmer basis. We find that the extrapo- 
lation of the step scaling function, computed with the 0(a)-improved action, to 
the contiuum need again quadratic terms in the fit. 

We refer to ref. U for more of the details about the calculation which we here 
only shortly summarize as follows. We calculate the renormalization constant of 
the twist-two non-singlet operator defined by: 

OfJ.» n = C^^h^ D, 2 ■ ■ • 3 M yV(z) + trace terms . (1) 

The basic ingredient for the reconstruction of the non-perturbative scale depen- 
dence of the renormalization constants of the above operator is the finite-size step 
scaling function o z defined by: 

Z{sL) = a z (f(L))Z(L) , (2) 

where L is the physical length that plays the role of the renormalization scale 
and Z the renormalization constant of the operator, which is defined by: 

R (fi) = Z(l/a/i)- 1 O barc (a/L) . (3) 



2 



Z is obtained from the SF matrix element of the operator on a finite volume L 3 T 
normalized to its tree level: 



O bare (a/L) = Z(L/a)O tree . (4) 

The framework of the Schrodinger Functional, which describes the quantum time 
evolution between two fixed classical gauge and fermion configurations defined at 
time t = and t = T, has been used extensively in the recent literature || |7j to 
calculate non-perturbative renormalization constants of local operators. Among 
the advantages of the method, we only quote the possibility of performing the 
computations at zero physical quark mass and of using non-local gauge invari- 
ant sources for the fermions without need of a gauge-fixing procedure. In our 
particular case, we exploit both features. Our observable is defined by M: 



f 2 (x = L/4) /ff 2 (x = L/4) 



Ti I \ Vf: 



(5) 



tree 



with f 2 given by 



/ 2 (* ) = -a^e^-^-^x)^ 5 2} V^(x)C(y)rVc(z)> (6) 



1 V, s ~ 1 , ,„1 



and fi by 



-3 ^3 



h = -a 12 E (C'(v)^c'(w)C(x)^C(y)) (7) 



y,z,v,w 



where ( = 5/5ip c and ( = —6/6ip c are the derivatives with respect to the two- 
component classical fermion fields and if) c , respectively) at the boundary x = 
0, while (' and (' are the corresponding derivatives at the boundary xq = T. The 
projection on the classical components is achieved by the projector P± defined by 
|(1 ±70). On the boundaries the theory possesses only a global gauge invariance 
that is preserved by the quantities defined above. The values of xq (set to T/4) 
and of the non-zero component of the momentum p x (set to 2ir/L) are both scaled 
in units of L, which therefore remains the only scale besides the lattice spacing a. 
The quantity f\ serves as a normalization factor that removes the wave function 
renormalization constant of the ( fields in order to isolate the running associated 
with the operator in eq. ([!]) only. 
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In order to compute the running of Z in eq. (|5]) its step scaling function oz = 
a z/ a h wr th 



has to be evaluated in the continuum. To this end, the physical size L (= T 
in our case) is kept fixed in physical units by keeping fixed the renormalized 
coupling constant in the SF scheme, renormalized at the physical volume scale. 
By increasing the number of lattice points while keeping the value of L constant, 
the lattice spacing in physical units is reduced and a continuum extrapolation of 
the step scaling function can be performed. 

The presence of a large non-zero momentum induces potentially large lattice 
artefacts in the numerator of the renormalization constant in eq. @. Indeed, 
in ref. for this reason, we had to extrapolate independently the step scaling 
functions o z and for cr^, with a quadratic and a linear fit in the lattice spacing, 
respectively. 

As already mentioned, the calculation done with a non-perturbatively improved 
action, without improving the boundary terms and the operator (whose matrix 
element in our case is tree-level-improved) does not eliminate all linear lattice 
artefacts and certainly does not eliminate the quadratic artefacts, if present with 
the Wilson action. Therefore, also with the clover action we have done indepen- 
dent extrapolations for <j z and for a/ 1: which are compared with our previous 
results using the Wilson action in tables 1 and 2 for the different values of L (or 
equivalently g 2 {L)) considered in ref. [[[]]. Within the errors, there is full con- 
sistency between the two sets of results obtained with different fermion actions, 
showing the universality of the continuum limit. Given the consistency of the 
extrapolated values of the step scaling functions from both fermion actions, we 
performed a combined fit, using both data sets obtained with the two fermion 
actions. In the combined fit we demanded that the extrapolations go to the same 
continuum value of the step scaling function. In figures [l] and |] we report results 
at fixed non-zero lattice spacing used for the extrapolation and the correspond- 
ing combined fit: while the linear dependence upon the lattice artefacts for er^ 
is almost eliminated, the quadratic terms for a 2 are still needed, as expected. 
In tables [I] and |2| we report the \ 2 P er degree of freedom (d.o.f.) values of the 




f 2 (x = L/A) 



tree 
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combined fit with those of the original fit of ref. [JTJ] and a fit for the case of the 
improved action alone. The x 2 values are totally comparable, as are our final 
values for the continuum step scaling function. Note that the error compared 
with the ones given in ref. |3J gets decreased by about a factor y/2. 

We have shown that the continuum extrapolation of the evolution of the renormal- 
ization constant corresponding to the average momentum of non-singlet parton 
densities shows the expected independence upon the lattice action used. Over the 
range of scales explored, the evolution is adequately described with a three loop 
expression for the anomalous dimensions JIJ. Further simulations that extend the 
running into the perturbative region, where a good description in terms of the 
known perturbative evolution can be made, are under way. 
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9 K L ) 


<j fa (Wilson) 


-NX 2 

X 


(j/^impr) 


-NX 2 
X 


<j f ± (combined) 


-NX 2 

X 


1.8811 


0.912(7) 


2.61 


0.910(7) 


0.06 


0.911(5) 


0.90 


2.1000 


0.897(9) 


1.17 


0.902(8) 


3.22 


0.900(6) 


1.54 


2.4484 


0.886(9) 


6.90 


0.884(7) 


0.17 


0.885(6) 


2.36 


2.7700 


0.866(10) 


0.52 


0.871(10) 


1.58 


0.869(7) 


0.73 


3.48 


0.837(11) 


0.39 


0.815(10) 


0.01 


0.825(7) 


0.86 



Table 1: The values of the step scaling function er^ are given extrapolated to 
the continuum. We compare the cases for Wilson fermions (2nd column), using 
the 0(a)-improved action (4th column) and from the combined fit (6th column). 
The running coupling is computed in the SF scheme. We give also the x 2 /d.o.f. 



9 2 (L) 


<7£ (Wilson) 


x 2 


o"^(impr) 


x 2 


U£ (combined) 


x 2 


1.8811 


0.917(8) 


0.19 


0.890(7) 


0.20 


0.902(5) 


2.55 


2.1000 


0.891(9) 


0.20 


0.891(9) 


0.89 


0.891(6) 


0.36 


2.4484 


0.891(9) 


1.01 


0.884(9) 


0.03 


0.887(6) 


0.43 


2.7700 


0.854(13) 


0.01 


0.858(13) 


0.002 


0.856(9) 


0.02 


3.48 


0.855(12) 


6.44 


0.840(13) 


0.04 


0.848(9) 


2.42 



Table 2: The values of the step scaling function are given extrapolated to 
the continuum. We compare the cases for Wilson fermions (2nd column), using 
the 0(a)-improved action (4th column) and from the combined fit (6th column). 
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Figure 1: Continuum extrapolation of Of x using a combined fit at all values of 
g 2 of table 1, with g 2 increasing from top to bottom. Circles and the dashed line 
correspond to using Wilson fermions. Squares and the solid line correspond to 
using the 0(a)-improved action. The fit for er^ is linear in both cases using the 
three points with smallest a/L. 
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Figure 2: Continuum extrapolation of 0% using a combined fit at all values of 
g 2 of table 2, with g 2 increasing from top to bottom. The fit for o % is quadratic 
using all data points. Notation is as in fig. 1. 
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